(7 " ARR Apr. 1942 




NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 



WARTIME REPORT 



ORIGINALLY ISSUED 

April 1942 as 
Advance Restricted Report 



A METHOD OF SHEAR-LAG ANALYSIS OF BOX BEAMS 
FOR AXIAL STRESSES, SHEAR STRESSES, AND SHEAR CENTER 



By Oscar Erlandsen, Jr., and Lawrence M. Mead, Jr. 
Grumman Aircraft Engineering Corporation 



NACA 



WASHINGTON 

NACA WARTIME REPORTS are reprints of papers originally issued to provide rapid distribution of 
advance research results to an authorized group requiring them for the war effort. They were pre- 
viously held under a security status but are now unclassified. Some of these reports were not tech- 
nically edited. All have been reproduced without change in order to expedite general distribution. 



W-33 



1176 01354 4649 
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EOR AXIAL STRESSES, SHEAR STRESSES, AND SHEAR CENTER 
By Oscar Erlandsen, Jr., and Lawrence M. Mead, Jr. 



SUMMARY 



A practical and relatively rapid method of compensat- 
ing for shear lag in box— beam analysis, with accuracy suf- 
ficient #for design purposes, is presented. Ef f ectiveness 
curves for box-beam elements are derived for an ideal, 
symmetrical structure. Application of the ideal curves to 
practical structures is described in detail using, as ex- 
amples, the center section of the wing for an airplane des- 
ignated as A and the unsymmetri cal D-beam of the wing 
for an airplane designated as 3. 

Tabular computation forms for rapid, accxirate calcula- 
tion of axial stresses, shear center, and shear stresses 
for a beam with shear lag are included. Results of analy- 
sis are compared with test stress distributions. Analyti- 
cal methods of checking the shear-lag curves by the use of 
the principle of consistent deformations are illustrated 
as a further indication of their reliability. 



INTRODUCTION 



Bending stresses in box beams with .thin cover sheets 
usually differ from those calculated by the ordinary bend- 
ing theory. This difference is caused primarily by shear 
deformation in the flanges or cover sheets. In order to 
compensate for differences from the theory, methods of 
shear-lag analysis have been developed. 

Kuhn (references 1 to 5) has done much mathematical 
and experimental work on idealized structxires under sim- 
ple loading conditions and has developed two methods of 
shear-lag analysis - the substitute single-stringer method 
and the shear-fault method - that approach the problem by 
means of arithmetic integration. This approach gives fair- 
ly accurate results for beams that are fairly symmetrical. 

^~.A S * ne beam becomes more unsymmetrical , however, with 
possibilities of changes of taper, camber, variable sec- 
tions, discontinuous stringers, holes, and indeterminate 
boundary conditions and as the loading pattern becomes less 



2 



simple, the NACA methods "become more difficult to apply. 
The computation is very tedious and of dubious accuracy 
where there are excessive approximations required of the 
designer in idealizing the structure so that it will fit 
the method. Furthermore, the NACA methods require two 
separate analyses for beam and chord "bending, which con- 
siderably lengthens the calculations. 

From design considerations, it is desirable to have a 
simple flexible method of shear-lag analysis that can be 
tailored to the individual conditions of any wing beam. 
It should require a minimum of idealization of sections so 
that irregular beams may be readily analyzed and should, if 
possible, enable a simultaneous consideration of beam and 
chord bending. It should be possible to check the results 
of the analysis in a simple straightforward fashion. 

The method of shear-lag analysis herein presented was 
developed by the Grumman Aircraft Engineering Corporation 
in the course of recent designs and, it is believed, incor- 
porates these desirable features. 

OUTLINE OF METHOD 



Essentially, the method consists in applying a correc- 
tion for shear lag to the areas of a beam cross section 
and then analyzing the substitute beam by the conventional 
bending theory. A check for the shear-lag correction can 
be made if desired. 

(l) Shear lag reduces the effectiveness of cover sheet and 
stringers for carrying bending stresses. Consequently, 
stresses in the cap strips will be raised and stresses in 
the stringers will be reduced from those computed by the 
ordinary beam theory. If the stringer stress at a section 
decreases 20 percent and the cap-strip stress increases 10 
percent from the theory, the effectiveness factor E' for 
the stringer will be 100 - 10 - 20 = 70 percent, referred 
to the beam theory cap-strip stress as 100 percent. 

Effectiveness curves are drawn in chordwise and span- 
wise directions for each stringer and wing section by tak- 
ing mathematically ideal curves and altering them arbitrar- 
ily to suit local conditions in the structure. Care must 
be taken that the location of the centroid of the effective 
material, and hence the shear center, varies in a smooth 
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spanwise curve. Common sense , cross-plotting, and a few 
rough checks will give surprisingly good results even with 
little experience in drawing such curves. 

(2) Stringer and effective sheet areas are reduced by the 
effectiveness factor E', and section properties are com- 
puted giving reduced "effective" values for I xx and I zz « 

(3) Axial stresses are computed "by use of the method pre- 
sented In reference 6 for unsymmetrical "bending, modified 
to take the effectiveness factors into account. These 
stresses should give the actual stringer and cap-strip 
stresses compensated for the shear-lag effect. 

(4) Shear stresses in the skin and webs that go with the 
axial stresses found in (3) can then "be accurately com- 
puted "by the new method. The shear center can also be ac- 
curately located and its position at the various wing sec- 
tions can be checked to see that it falls on a smooth 
spanwise curve. 

(5) The accuracy of the original selection of E' can be 
proved by utilizing the principle of consistent deforma- 
tions. Any two points on the wing surface can be selected 
and the deformation of one of them from the other can be 
computed by two or more different paths, using the axial 
and th.9 shear stresses obtained in the analysis. These de- 
formations should agree if the values of E* were reason- 
ably correct. 

COMPARISON WITH NACA METHODS 



The main difference between the effectiveness-factor 
method and the NACA methods of shear-lag analysis is the 
point at which engineering .judgment is required. . The NACA 
calls for judgment in idealizing the section before shear 
lag is considered. The effectiveness-factor method re- 
quires an estimate of the shear-lag effect at the outset. 
It is believed that in the design of new wings, the 
effectiveness— fact or method is more direct, more rapidly 
applied, and comparably accurate. 

In the NACA methods it is easy for the average person 
to lose sight of the physical significance of the problem 
amidst the maze of coefficients and unfamiliar formulas 
and to make time-consuming mistakes. The effectiveness- 
factor method obviates this possibility, it is felt, by its 
simplicity and flexibility. 
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(1) Once the effectiveness curves have "been draxvn, the 
analysis follows the standard beam theory. Both "beam and 
chord "bending are considered simultaneously, as described 
later. 

(2) The method is perfectly applicable to any shape of 
"beam "because the areas of the existing structure are 
changed rather than the shape of the "beam "being idealized 
without changing the areas, as is done in the HACA method. 
It is believed that changing the area is an easier task, 
particularly for complex uns,ymmetrical sections such as D- 
beams, where idealization of shape would be very difficult. 

(3) Local irregularities and changes in section can be 
readily considered by suitably drawing the effectiveness 
curves. Material outside the main beam may also be con- 
sidered in the calculations if it is believed that it 
carries bending stress. This consideration would be dif- 
ficult by the NACA method. 

(4) Effectiveness curves give a handy tool with which to 
design stringer skin-rivet connections. 

Analyses by both the NACA' and the effectiveness- 
factor methods vrere made for a wing section of airplane A 
and were compared with test results. Remarkably good 
agreement was shown for both methods, but the effectiveness- 
factor method showed more adaptability to local conditions 
and required less computation. 



EFFECTIVENESS -FACT OR CURVES 



The nature of effectiveness factors and their appli- 
cation to shear-lag analysis will now be developed and 
illustrated for a simple case. A mathematical derivation 
of formulas for the values and the shape of effectiveness- 
factor curves in both spanwise and chordwise directions 
will be presented for an idealized beam. The method of 
considering beam and chord bending simultaneously with 
shear lag included will be outlined, and the relation be- 
tween effectiveness curves and the location of the shear 
center in un symme t r i cal ' sect i on s will be discussed. 
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Nature of Shear- Lag 

Shear lag and how to compensate for it can be illus- 
trated by consideration of figure 1, which represents the 
bottom surface of a wing beam from the plane of symmetry 
of the airplane to the hinge fittings for the outer panel. 
Here two loads P act on the cap strips. 

The usual theory assumes that the panels will elongate 
as shown in figure l(a), which would require uniform stress 
across the panel at all .sections, even across section D. 
Thi.s assumption is obviously impossible beca-use the center 
stringer can have no stress or axial load at its outer 
free end. 

Actually the panel will deform as in figure l(b) ow- 
ing to the shear deflection or "lag" in the sheet. The 
thinner the sheet the more the lag. At the free end the 
stress in the stringer will be zero. As one progresses 
along the stringer toward the root, the load P starts to 
spread out into the sheet in the form of shear stresses 
that load up the stringer. 

At the root for an infinitely long panel the stress 
in the stringer would be equal to that in the cap strips 
and there would be no shear stress in the sheet and no 
shear lag. 

It is important to realize that the stringer can be 
loaded only by shear from the sheet. If the panel is 
short or the sheet very thin, the sheet will not be able 
to transmit enough load to the stringer to stress it at 
the root to the same amount as the cap strip (as would be 
figured from the ordinary beam theory). This case is usu- 
ally true for wing box beams. 

Hence, it is seen that the stress in the stringer 
varies from zero to some maximum value at the root; this 
value depends on the dimensions of the beam. Compared 
with the uniform stress of P/A calculated from the beam 
theory, the stringer is operating at a reduced stress, or 
reduced effectiveness, that varies on some smooth curve 
from 2' = 0 at the free end to E 1 = maximum at the 
root. Likewise the cap-strip stress will vary from a max- 
imum value of the free end to a minimum value of the root. 

These deviations from the beam-theory stress may be 
obtained from figure 2. 
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At section 0 of figure 1(b), let it be assumed that 
with shear lag the stringer stress will "be half of that 
computed "by the P/A formula (fig. 2). The value of E» 
for the stringer would he 0,50 and for the cap strips 
would be slightly greater than 1.00. It is convenient, 
however, always to assume the E for the cap strips to be 
unity. . 

The effectiveness E' for a stringer will be defined 
as the effectiveness of the stringer for carrying stress 
as compared with the cap strip. This definition will make 
E • for the stringer slightly less than 0.50, say 0.43. 
If the effectiveness factor for each member is applied to 
the area of the member, an effective area is obtained. 

Effective area = (3x1 . 00) +( 2x 0 .43 ) +(3xl .00) =6 .86 square inches 

p 

Shear-lag stress in any member = X E* 

^effect 

20000 

For cap strip " 6 3 g x 1.00 = 2920 pounds per square inch 



For stringer 2 Q°?2 x 0.43 = 1254 pounds per square inch 

6.86 

Beam-theory stress = P/A = 20 9°° = 2500 pounds per square inch 
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Shear-lag stringer stress 1254 

— = = 0.50 

Beam-theory stringer stress 2500 

This result checks the original assumption. It is pos- 
sible, therefore, to take shear lag into account by select- 
ing effectiveness values for each element of the cross sec- 
tion. The component areas are corrected by these factors, 
and from there on the calculation is similar to that for 
the ordinary beam theory. The actual value of the effec- 
tiveness factors selected will be governed by the mathemat- 
ical derivations to be discussed later. 



Mathematical Derivation of Spanwise Curves 

A box— beam cover is statically indeterminate internal- 
ly to the first degree; therefore, with two static rela- 
tionships and one elastic relationship, it can be solved. 
The symbols are defined in appendix A. From these three 



fundamental equations, differential equations for f p , fj,, 
and f g can be derived in terms of the physical constants 
and the loads of the problem* These differential equations 
can be solved mathematically for simple loadings and sym- 
metrical sections. The ratio of fj, to fp will be the 
effectiveness factor for the stringer member in terms of 
the cap strip. 

The sketch of figure 3 shows the idealized beam with 
applied moments and shears. The material in the bottom 
of the beam can be considered to be concentrated entirely 
at the bottom flanges. The origin of the coordinates is 
taken at the tip. (See reference 1.) 

Ef f ectiveness.-factor curves for the idealized beam 
sketched in figure 3 will be mathematically derived. The 
following assumptions are made in the derivation: 

1. Rib,s, and hence beam, have infinite chordwise stiff- 

ness so that chordwise strains may be neglected 

2. Diagonal tension in the sheet is taken care of in 

factor Gr 

3. Sections and loads are symmetrical about the yz plane 

Fundamental equations .- Consider the static equilib- 
rium of the element dy (fig. 4 and reference l). It will 
be seen that 

dP L = dS L (la) 



= _ S 4^ « aSfl • db) 



dP;E , = _JL d s 0 

From assumption (l) and figure 5, it follows that the shear 
stress 



f_ = £ (AF - AL) 
s b 



but 



■I 



. £.-■ 

y=l 



and 
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fp dy 



i 

°y=L 

and, "by differentiation, 

df s = - ^ (f F - f L) ^ (lc) 

Effectiveness factor .- Equations (la), (it)), and 
(lc) can "be combined into four differential equations for 
f™, f L , and f g . (See appendix B for detailed derivations.) 

For the beam in figure 3 the equations for fp and f^ 
can he solved to he 

f,,JL (i3 cos& K (L - r> 1 ( 2a ) 

Afph L Ap cosh :CL J 

fi - cosh K (L - y) 1 (2b s 
L cosh XL J 



*L = 



M 

Arph 



where II/Arjtla is the ordinary beam-theory stress and the 
term in the bracket is the effectiveness factor of each 
referred to the beam-theory stress as unity. Since the 
E 1 for the' flange, or cap strip is taken as 1*00 for con- 
venience, the actual E' for the stringer referred to the 
cap strip will be 

1 _ cosh K (L-y) 

E« = It = " c °s h KL (3) 

1 + ^k. cosh K (L-y) 
Ay cosh KL 



where K, the shear-lag parameter, equals 

and the other symbols are defined in appendix A. Equation 
(3) is plotted in figure S for various values of K and L 

Limitations and Application of Spanwise Curves 
Mult i stringer beams .- Equations (2b) and (3) are for 



/9± C_L + J.) 

J Eb \kv At, / 



a single-stringer "beam. When there is more than one 
stringer, the usual case, the stringers "between the cap 
strip and the shear center should "be lumped into a single 
stringer at the area centr old of the stringers, a distance 
b . .fr.o.m the. .cap. strip., -bef ore . K- is. -computed. -Ithi-s meth- 
od will give a good approximation of the average E ' for 
the stringers for a fairly symmetrical box beam. The chord- 
wise distribution between these stringers will follow the 
shape of the chordwise curves discussed later. 

For very unsymmetrical "beams, such as D-beams, values 
from equations (2b) and (3) will be only very approximate, 
but the general shape of the curves will always hold and 
the values can best be obtained from consideration of the 
location of the shear center and the center of gravity of 
the effective material. 

Variable sections . - Equation (3) is strictly true 
only if, at each spanwise station, the following condi- 
tions hold: ; 



(a) Constant K ; 




(b) Constant Al/Aji and no camber 

(c) Constant beam-theory axial stress M/Ajh ■: 

(d) Constant vertical-web shear stress Sw/ht 

Conditions ' (c) and (d) are sought for in practical . 
wing design and condition (b) is usually fairly nearly 
true. Condition (a) is seldom satisfied 'for a tapered 
variable-section wing but, by the use of an average value 
for the panel, the error will not be large and judicious 
fairing of curves will give results accurate enough for 
design purposes. '• • 

Riveted .joints .-. Tor a constant-section stringer, 
the slope of an effectiveness curve such as shown in fig- 
ure 6 represents to some scale the change in axial load 
per inch of span. This value is the design load for rivet 
spacing between stringer and skin and wiTl be critical at 
the free-end of the stringer where the slope is steepest. 
Since for most wings the curves will not be precise, the 
1.20 fitting factor should not be omitted in the design of 
rivets based on effectiveness curves. The curves of fig- 
ure 6 are critical for stringers but are not critical for 
rivets; the 1.20 fitting factor will cover any divergence. 
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Mathematical Curves for Chordwise Effectiveness 

In references 1 and 2 analytical solutions have "been 
made that show the chordwise distribution for ideal beams 
with a cover sheet without stiffeners (or infinitely close- 
ly spaced stiffeners). This distribution follow* a hyper- 
bolic cosine law 



c o sh 



(»!) 



* cosh Yb - 

where fp is the stress in the cap strip (x = b) and 1 
is an auxiliary parameter obtained from the relation 



tanh Yb 



=4^ (5) 



Yb f F 

where ^i/^s 1 * s the ratio of the stress in the stringer 
material located at the stringer force centroid to the 
stress in the flange. Figure 13 of reference 1 is a curve 
presented to facilitate solution of equation (5). The 
ratio of f/fj, from equation (4) will be the effectiveness 

of any stringer in terms of the cap strip or flange. Curve 
of f/fj, for various values of Yb or f ^/f p are pre- 
sented in figure 7. 



Limitations and Application of Chordwise Curves 

In general, the mathematical curves of chordwise effec 
tiveness are not strictly accurate, because 

(1) Most wings have a finite number of stringers and will 
not follow equations (4) and (5) exactly 

(2) Local conditions, dissymmetry of section, etc., will 
vary curves 

The . theoretical curves, however, will give the general 
shape of chordwise curves and also will give a point of de- 
parture from which to correct for the local conditions of 
the problem. 
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Curves for Simultaneous Beam and Chord Bending 

The previous discussion has dealt with effectiveness 
curves for bending in only the "beam direction. Wing beams 
a're subjected to both beam and chord bending. For the box 
beam shown in figure 8 the chordwi se-ef f ectiveness curves 
for beam and chord bending are shown and it will be noted 
that they are not the same. For a precise analysis, this 
difference means that there are two different beams with 
different areas corresponding to the two effectiveness 
curves. The beams must be separately analyzed for beam 
and chord bending, which would be very tedious. 

If a single value of effectiveness for both beam and 
chord bending can be selected, the beam may be analyzed by 
the unsymmetrical-bending formula. This will be done by 
using the beam-effectiveness curve as the true total effec- 
tiveness of the beam. The error introduced will be small 
for the following reasons: 

(1) The chord moment, at the most, is only about one- 
third of the beam moment and it is usually less than this. 

(2) For the beam of figure 8, where the beam and the chord 
effectiveness of members b and c near the chord neutral 
axis differ markedly, the stresses due to chord bending 
will be small and the errors even smaller. This result is 
true in most box-beam sections. 

(3) Material back of the main beam is usually ineffective 
in beam bending but quite effective in chord bending. In 
order to take care of this condition, the effectiveness 
factor E' will be selected from consideration of chord 
bending and the z coordinate of the material will be as- 
sumed to be zero, so that the terms will drop out of the 
beam-bending calculations just as though the beam effec- 
tiveness were zero, but Will not drop out of the chord cal- 
culat i ons . 

(4) Where beam effectiveness is very low and chord effec- 
tiveness is high, the value of the effectiveness factor 
B' can be slightly increased over the value for beam ef- 
fectiveness in order to compensate in part. 

A little thought and cart, i n . e valuat i ng the relative chord 
and beam effects will give effectiveness' curves that are 
sufficiently accurate for design purposes. 
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RELATION OF SHEAR CENTER TO EFFECTIVENESS CURVES 



In symmetrical sections the ideal curves of figures 
6 and 7 are usually sufficient guide for obtaining fairly 
accurate results. In unsymmetri cal sections guch as D- 
beams, however, the ideal curves give a guide as to the 
shape of the curves hut are of little service in selecting 
the values, because it is nearly impossible to get an accu- 
rate value of K with which to enter the curves. The lo- 
cation of the shear center and the centroid of the material 
effective in bending are used as guides in drawing effec- 
tiveness curves for such beams. 

Shear-center location .- The shear center, s.c, of a 
section shall be defined as the point at which any load may 
be applied to the section without tending to twist the sec- 
tion. Only pure bending will occur. (It should be noted 
that the sheer center is not necessarily the same as the 
elastic axis about which the section twists.) For a sym- 
metrical section the shear center will coincide with the 
centroid of the section. For unsymmetrical sections, how- 
ever, the resultant of the shear forces on the section 
does not pass through the centroid. (See fig. 9.) 

The point s.c. is so located that the moment of the 
shear forces about s.c, is zero. Any load applied 
through point s.c. will cause pure bending in the direc- 
tion of the load. Reference 6 shows that for. a D-beam 
section the shear center s.c. lies as shown in figure 9(c), 
inside the closed torque box and not far from the centroid. 
The heavier the vortical web becomes with respect to the 
skin sheet thickness, the closer the shear center will move 
toward the web. 

Shear center with shear lag .- The presence of shear 
lag is taken care of by a reduced effectiveness for the 
area of material other than at the cap strip. In a D-beam 
this reduction will move the centroid of the effective area 
closer to the web than was the centroid of the total area 
and hence closer to the shear center. For sections with 
closed torque boxes it will be assumed, therefore, that the 
shear center and the centroid of the effective material 
coincide and fall within the torque box. This assumption 
is only approximate but the error introduced is small. The 
assumption is a useful one, because rules that apply to the 
location of the shear center can be applied to the centroid 
of the effective areas and will serve as a control on the 
shapes of the. effectiveness curves selected. 
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This law is used "by drawing effectiveness curves such 
that the center of gravity of the ef f ective -material; -_f<Xl- c I 
lows a smooth easy spanwise curve such as the shear center 
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give fairly definite values for the effectiveness at any 
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Sample effectiveness curves for wings of two '"airplanes , 
airplane A and airplane B will he drawn. Calculations for 
the corresponding stresses were made f.or^-th-e. )v/xng-of ^air- 
plane A (not included herein) and comparison with test re- 
sults will he made. For the wing of airplane B, sample 
calculations will be included,-. t q illustrate the, standard 
forms th a t~have been "developed for rapid" wing analysis? 

r- i , Wing Beam . of Airplane A < _, : r ^v.,.;>o ft o ? 

Center section .- ,.5hear 7 lag effectiveness, curves- ^ill,- 
be drawn for the wing .of ..airplane A.. from stations .18 to 110. 
Figure 11 shox^s a typical cross section and plan; Inboard of 
station 13 the section will be consider e$ _stiff ? enough to 
neglect shear lag, becausTe-Vf the 0.081- inch doubler plate 
(later altered to 0.102 inch but here'kept at 0.081 inch 
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for comparison, with tests)* Outboard of station 110 there 
is a large hole in the top that will he considered. Mate- 
rial ahead of and hack of the main "box "beam will he consid- 
ered ineffective and will he neglected. Effectiveness fac- 
tor 15 1 is plotted for the top cover in figure 12. 

The effectiveness in chord "bending would he unity 
throughout the c-over hut, as stated previously, the stress- 
es due to chord heading will he small and the curves will 
he drawn according to "beam "bending only. Because of symme- 
try, the shear center of the original section will be ap- 
proximately at the center of gravity of the section. Shear 
lag will not appreciably move the shear center and a de- 
tailed study of these motions will be unnecessary in this 
simple structure. (See airplane B D-beam treated later 
for the case in which shear-center location is significant.) 

Effectiveness, factor determination for ton surface .- 
The effectiveness factor for the top surface can be deter- 
mined as follows: 

1. Stringer D - 

Find spanwise-eff ectiveness distribution for stringer D 

At station 18, E' =>0 

because of hole beyond these 
At station 110, E ' = 0 points 

To approximate E 1 at station 66, compute for this 
section 

v-s _ Gt (l_ . 1 ^ 
Eb \L P A t ' 

~ = 0.40 t = 0.156 b = 8.00 Aj. = 1.004 A L = 8.702 

These values were obtained from idealization of sec- 
tion 30-94 (fig. ll). Consider one-half width of beam and 
take b as the distance from front beam to center of grav- 
ity of stringers and cover plate. Ap = elements (l) + (2) 
+ i(3), A L = elements £(3) + (4) + (5) + (6) + (7) + £(8) 

K 2 = 0» 4 Q x 0*156 x 1 + 1 

8.00 1.004 2.702 

= 0.0078(0.995 + 0.370) = 0.01066 
K = 0.1032 



Prom the curves of figure 6 the ef f ectiveness at oen- 
ter of span is 0.987 for K = 0.10. Therefore, assume that 
E 1 = 92 percent at station 6S for stringer D. This value 
is considerably lower than 0.99 hut is justified because: 

a. The quantity K is very approximate owing to irreg- 
ularities of section, camber, etc. 

b. Curves of figure 6 are for a single stringer at cen- 
ter of gravity of stringer material on one-half of center 
line. For stringer at center line of total section, B' 
will be less, following smooth curve of figure 7„ 

Inboard of station 66, curve D can be assumed to fol- 
low very generally the shape of the idealized curve and is 
held high in order to take care of the increased stiffness 
from stations 18 to 30 because of the doubling plate. Out- 
board of station 66, the curve has the idealized shape but 
is lower than the inboard curve because of the thinner 
cover sheet. 

2. Stringers B and E - 

Stringers B and P can be considered fully effective 
over almost the entire length of the beam, because of 
their proximity to the flanges and their continuity through 
the plane of symmetry of the airplane. At station 126 the 
effectiveness will be zero because the hinge fittings of 
the outer panel are attached to the beams only. The effec- 
tiveness curve will stay very high until quite near sta- 
tion 12S, so that in the region from stations 94 and 110, 
E' will drop only slightly. 

S. Ohordwise Plots - 

With tentative E' curves plotted spanwise for string- 
ers A, B, and D, chordwise plots through these stringers 
may be made using the shapes of ideal curves and with fig- 
ure 7 as a guide. 

4. Stringer 0 — 

Erom the chordwise plots of step 4, the spanwise 
curve for stringer 0 may be obtained. Hote that this curve 
Is similar to that of stringer D but is slightly higher. 
The shape of this curve is a cross check on the accuracy 
of the chordwise plots. 
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5. Stringer E - . 

At station 18, £' will "be very nearly equal to 1,00 
"because of continuity. At station 110,. £' equals 0. 
Therefore, the curve will he as shown in figure 12. At 
the outboard end it will coincide with the curve for 0 
because of the symmetry. The chordwise plots may now he 
completed. Bottom-surface effectiveness curves were also 
drawn but are not shown. 

Comparison of analysis with test re . sults .— Final 
stress values ohtained hy the effectiveness-factor method 
of analysis from the curves of figure 12 are plotted in 
figure 13 as a chordwise plot and are compared with test 
readings taken with Huggenberger strain gages on the ac- 
tual wing surface. These test results are as reliable as 
test results can be expected to be, for there was no buck- 
ling in the skin nor undue distortions during the observa- 
tions. Curves for the stress by the beam theory and by 
the NACA suhstitute single— stringer method of shear-lag 
analysis are also presented for comparison. 

Agreement between test values and the ef f ectiveness- 
facter method is generally very good, and it is easily ap- 
parent how to alter the effectiveness curves to improve 
the agreement. It was assumed that stringers "B and T 
were fully effective and stringer E nearly so. This as- 
sumption was conservative. If the values of E T of these 
members were dropped slightly below 1.00, the agreement 
would improve. The curve as drawn, however, is sufficient- 
ly good for design purposes. 

Note that, with the test results, compressive stresses 
tend to be lower at the front beam and higher at the rear 
beam than computed "because of the effectiveness .of the ma- 
terial outside the box beam, particularly for chord bending. 

The HACA method for just beam bending indicates that 
good agreement could be had if the chord effect were super- 
posed on it. The necessity for idealizing the section ac- 
counts for the symmetry of the curve and the lack of local 
irregularities that are picked up by the effectiveness- 
factor method. The ordinary beam theory shows good agree- 
ment only at the front beam and emphasizes the need for a 
quick and relatively accurate shear-lag method. The 
effectiveness-factor method seems to fulfill this require- 
ment admirably. 
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Ving D-Beam of Airplane B 

The D-beam in the wing of airplane B is quite unsymmet- 
rical and hence compensation for shear lag is more compli- 
cated than for the ordinary "box "beam. In addition to hav- 
ing smooth effectiveness curves in chord and span directions, 
based on the idsal curves of figures 6 and 7, the location 
of the shear center and the centroid of the effective mate- 
rial must he made to follow a smooth curve, 

A complete shear-lag analysis of this wing has been 
made hy the Grumman Aircraft Engineering Corporation by the 
use of the effectiveness-factor method. In the present re- 
port a detailed description of the reasoning "behind the ef- 
fectiveness curves is given. Certain modifications of the 
curves are suggested, and new calculation forms are pre- 
sented to speed up future wing analyses. 

In the analysis separate curves for "beam and chord ef- 
fectiveness E' B and were drawn and used in the 
unsymmetrical-bending formula. This method is not strictly 
accurate mathematically and instead the method for a single 
£' is recommended. In the case of the wing for airplane B, 
the error is insignificant "because the product of inertia 
is very small at most sections, which would not necessarily 
"be true of all wings. 

Location of shear center .— The "best place to start 
the shear-lag analysis is at the hinge fittings "between 
the center section and the outer panel of the wing where 
the location of the shear center and loads are completely 
determined. (See figs. 14, 15, and 16.) For positive 
lift, all the tension load is carried at C (fig. 17); 
therefore, the shear center lies on a vertical line through 
C and the compression load is divided "between A and B 
inversely as their distance from the shear center line. 

Inboard of the hinge fittings, the material in the 
"bottom cover forward of the cap strip C will start to 
take some of the compression load as the sheet transmits 
shear around from A and forward from C. This action 
moves the center of the compression load forward and with 
it the shear center. At the plane of symmetry of the air- 
plane the shear center will "be the farthest forward but, 
"because the shear-center line must vary in a very easy 
smooth curve, It will not come as far forward as though 
there were infinite length to the center section. 
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Outboard of the hinge fittings the presence of the 
wheel well in the "bottom surface prevents the compression 
load from distributing forward. This result keeps the 
shear center in line with the lower cap strip until out- 
board of the wheel well. Prom there to the tip, the shear 
center will move forward as the compression loads spread 
forward, until at the tip the shear center will coincide 
almost directly with the location according to the beam 
theory. The motions of the shear center are shown in fig- 
ure 14. A method of checking the location of the shear 
center will later "be described. 

Effectiveness-factor determination for ton surface.— 
If the assumption that the shear center and the centroid 
of the effective areas coincide is used as a guide and the 
curves are faired to look like the ideal as much as is 
practicable, the effectiveness curves in beam bending may 
be drawn. (See fig. 15.) 

She cap strip B will be considered to have unit ef- 
fectiveness throughout. 

1. Inboard of hinge fittings - 

The value of B 1 for stringer 5 will fall off rap- 
idly as the leading edge H and stringer 0 build up, 
in order to keep the centroid of the effective material 
back with the shear center . At the plane of symmetry 
(station 0) stringers C and H will have an effective- 
ness of about 55 percent because it is impossible to load 
them any more than by means of the shear in ^he sheet . 
Stringer G will come up a. little from a low point at 
station 44 as the shear center moves slowly forward. 

The ideal curves of figures 6 and 7 are tised to fair 
in the shape of the stringer curves. It is believed that 
the curve for stringer B would be improved if it were 
made to drop off to zero between stations 29 and 44, since 
it is impossible to stress the end at station 29 with more 
than the sheet can transmit in compression. 

3. Outboard of hinge fittings - 

The effectiveness of stringer G- will have to fall 
off very rapidly through the wheel well as stringers 0, 
D, E, F, and H build up from zero in order to keep the 
shear center and the centroid of effective area loads over 
the lower cap strip. Stringers C, D, E, and P will take 
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more or less ideal shapes, dying out to zero at each end. 
The maximum height of these curves will he governed by the 
condition that the shear center cannot move forward too 
fast outboard of the wheel well and by the rate at which 
the top cover sheet can transmit shear as shown by. the 
slope of the chordwi se-ef f ectiveness curves. 

Stringer H is kept very low until outside the wheel 
well because of the discontinuity of the torque box through 
this region and the nature of the curves for the under sur- 
face. Outboard of the wheel well H will gradually rise 
to unitjr near the tip as the shear forces in the sheet load 
it. Stringer G will drop to zero at station 277. 

3. Material back of cap strip - 

To neglect completely material back of the cap strip 
as ineffective in beam bending would introduce appreciable 
error, particularly for stringers near the cap strip, where 
it is necessary to transmit considerable shear in the sheet 
and riveted connections. In the analysis of effectiveness 
factors for beam and chord bending, the curves were as 
shown in figure 15» If a single factor E' is used, as is 
recommended for future curves, the E 1 for member A would 
follow the E' B of figure 15. The E« for the trailing 
edge material (element V) would be taken equal to E'q of 
figure 15, and the material would be assumed concentrated 
at the beam— bending neutral axis ( z = 0) so that it would 
enter the chord-bending "but not the beam-bending calcula- 
tions • 

4. Pairing in curves - 

The final curves are obtained by working between the 
chordwise curves and the spanwise curves, with the ideal 
shapes and the relative effects of beam and chord bending 
in mind. They should show the general behavior of the top 
sheet without consideration for the minor irregularities 
that must be taken care of with more detailed investiga- 
tions . 

Effectiveness-factor determination for bottom sur- 
face . — The lower cap strip will be considered 100 percent 
effective throughout in computing the effectiveness fac- 

• 4or -£or,.-tie bott om .surf ac^e,,^. 

1. Inboard of hinge fittings - 



20 



Stringer J will "be zero at the hinge station "because 
the fitting is only on the top chord. It will build up as 
shown in order to keep shear center hack where it should he. 

Spanwise plots I» and L 1 need no comment. In the 
chordwise plots of figure 16, however, that for station 44 
seems in error. The dotted line should he followed, not 
the solid one, since the hole outboard of station 44 pre- 
vents stressing the sheet at this point, even in tension. 

2. Outboard of hinge fittings - 

In the wheel-well region, stringer J is subject to 
local deformations caused by the discontinuity of the 
torque box through the wheel well. The inclusion of these 
effects in the general analysis \fould be very difficult; 
the effectiveness for member J was therefore taken at 
zero from the hinge to station 157 for the bending-section 
properties and the local effects were later superimposed 
there on . 

Springer L has very nearly the ideal effectiveness 
shape, but the ideal curves would indicate that it should 
not cone up from station 142 quite as steeply as is shown. 

3. Material back of cap strip - 

Effectiveness curves for the material back of the cap 
strip are drawn in a manner similar to that for the top 
surface . 

COMPUTATION J GEMS FOR BENDING STRESSES 



Tabular forms have been prepared to facilitate compu- 
tation of axial stresses and loads in wing beams; samples 
of the forms are included as tables 1 and 2. Table 1 
shows the calculation of section properties for the effec- 
tive areas. Columns 1 to 17 and computation of I X x» *zz» 
and I xz need no explanation. The symbols C lf C s , C 3 , 
and C 4 represent terms in the unsymmetri cal bending formu- 
ula altered to suit the effectiveness-factor method of 
analy si s . 

m V_^cixx E , x+ M 0 K - ,,. . Kl E«x + K 3 E'z 

3 I I - I II - I 

■ L XX- L ZZ X XZ J- XX. ZB XB 



The values of K x and K 3 in columns 26 and 27 for 
each flight condition are used in table 2 to compute the - 
"bending stresses and the axial loads in each member. 
These load values are used in obtaining the shear stresses 
in the section and need be calculated only for those condi- 
tions in which the shear stresses may be critical. 



SHEAS STRESS AND SHEAR CENTER 



As a rough check on the accuracy of the effectiveness 
curves, it would be desirable to be able to find the loca- 
tion of the shear center in a fairly rapid, straightfor- 
ward fashion. This shear center should correspond to the 
axial stresses and loads that were obtained from the curves 

The shear center may be found by computing the shear 
stresses in the sheet and the webs due only to bending and 
then by finding the point in the section around which the 
moments of these shear stresses are zero. It is an estab- 
lished principle that, when the axial bending stresses and 
loads in the members have been selected (whether rightly or 
wrongly makes no difference), there is only one combination 
of shear stresses that will maintain equilibrium. This 
principle has been discussed fully in references 6 to 10 
and will now be discussed in detail. The principle states 
that, for any set of axial loads, it is possible to obtain 
a corresponding shear center. The spanwise curve of these 
shear centers will then serve as a check on the effective- 
ness curves, and the shear stresses obtained can be used in 
the design of the webs, sheet, and riveted connections. 



Fundamental Equation for Shear Due to Bending 

Consider the equilibrium of the element shown in fig- 
ure 18 representing a portion of a box beam. Assume that 
the effective area of the sheet in bending is combined 
with the stringer. If q n is the running shear in pounds 
per inch at station n, then 

dP J 

where dP/dy is the average change in axial load per inch 
in the stringer over the elemental distance dy, and q x 
is the average shear per inch over section 1. 
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If there are several other stringers between sections 
1 and n (fig. 15), then the shear per inch at section n 



will he 



n 



«n - f £ + fa <?> 



or, in terms of shear stress where q = f e t , 



n « 

The term 2 ^ is easily obtained from the computed ax- 
ial loads. If a load curve is drawn for each stringer and 
cap strip, the slope of this curve at any station will be 
dP/dy. A less accurate method that can be used if desired 
is described in reference 7; it gives the average shear 
over a finite length of wing. Table 3 is applicable to ei- 
ther method. 

Derivation of Equation. for Shear at Cut Section, q x 

The final value oT shear flow or running shear q n 
at any point due to bending alone must be such that the 
total shear deformation around the section from section 
n back to section n. will be zero. Only one value of 
q,x will make q. n fulfill this requirement. From the 
method given in references 6, 8, and 9, q x may be ob- 
tained by the method of least work. 

Reference 10 presents a method of finding the total 
shear in a section, both bending and torsional t by a stat- 
ically determinate method without the least-work solution. 
Reference 7 has amplified the method for general applica- 
tion in a simple straightforward fashion that completely 
eliminates computation of the shear center. If the shear 
stresses were the only consideration, this method would be 
the most practical. As a check on the shear-lag effective- 
ness curves, however, it is here desirable to be able also 
to find the shear center in a rapid fashion, which can be 
done only by making use of the least-work solution and ob- 
taining the shear stresses, Just due to bending, without 
the superimposed torsion. The computation is not much 
longer than that shown in reference 7 and the additional 
information obtained is worth the difference. 



The least-work solution requires that 
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3f B 



= o 



where V is the total work due to shear deformation in the 
sect i on, 

integrated around the complete section when f 8 t ** q 8 , and 



But, from equation (7a) 



so that equation (8) "becomes 

L 



" \f f =n *» dS " 0 



0 

For any sectipn, G can "be considered a constant and 
will "be assumed to include the effect of diagonal tension. 
Therefore 

L 



/ 



f_ t, ds m 0 (9) 
8 n 1 



If equation (7a) is substituted in equation (9), 

. -X^i wY-£tkl:L.. 



or 
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q x = (10) 



0 *a 

If aquation (10) is substituted in equation (7), 




(11) 



Equation (ll) 1b perfectly general and is applicable to 
any single cell bean. Hote that q n is just the shear 
due to bending and assumes the shear loads to be acting 
through the shear center. Shear q n is also the average 
shear flow p'er inch over the spanvise distance dy. 



Double Oells and Other Shapes of Beams 

If the box girder has more than one cell or has mate- 
rial extending out beyond the single cell, the shears may 
be found by extending the same principle. A detailed der- 
ivation for finding the redundant shears in a double tor- 
sion box is presented in appendix 0. 

In the type of. section shovn in figure 20 the shear 
in the extended flanges is scatioally determined because 
at the free edge the shear is known to be sero and the 
shear at any section of the extension may be obtained by 
the relation 

q - S (12) 



The shear from equation (12) can then be directly added 
to the shears in the box, and the problem is reduced to 
that of a single torque box. 
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Location of Shear Center 

Once the shear stresses due to axial loads have "been 
determined the shear center may be located. If the ex- 
ternal shear loads are considered to act through- the shear 
center, the moment of the shear stresses about any point 
such as the center of gravity must be equal to the moment 
of the shear loads about that point, for static equilib- 
rium. 

The moment of the shear flow q about the center 
of gravity is equal to ,q times twice the cross-hatched 
area in figure 21. 



2dA = x n z n+i 



" x n+i z n 



For a section of large curvature, such as in the nose 
D, the area between chord and curve may be estimated, dou- 
bled, and added to formula (13), if increased accuracy is 
desired. 

Now resolve V and V into a single shear V and 



take moments 



(V X S + Y_ S ) d = Vd = Zq 2dA 



d = 



Eq 2dA 



(14) 



where d is the distance from the shear center to the 
center of gravity. Therefore, 



(15) 



v x _ v x S °. 2dA V x 2( l 2dA 
z s.c. - y V 3 = V 2 * + Y z ~ 



The location of the shear center for each station should 
plot in a smooth spanwise curve if the original effective- 
ness curves are fairly accurate. 
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Shear Due to Torsion 

Since Vj. and V z are not located at shear center, 
there will "be a net torque on the section about the sheai 
center. 

T s.c. = T 0 + r z (x c.g. - «..o. ) - V x (z c.g. + ^.cj 

The shear flow on any section due to this torque will 

he 

T 

q t = 2k 

This shear flow can then he added to that due to the beam 
and the chord tending to give the total shear flow and 
hence the shear stress at any point. 



Effect of Taper 

Taper in hoth depth and width will affect the shear- 
stress distribution. Taper in width can usually he safe- 
ly neglected. Taper in depth, however, is important in 
the design of the vertical webs, since the cap strips will 
carry a fair share of the applied shear load and relieve 
the vertical webs. 

Consider the tapered panel in figure 22. Take moments' 
about A and assume dy so small that V may be consid- 
ered constant over its length. 

T dy + Ph - dPh - Ph - Pdh = 0 

7 dy - dPh ^ Pdh = 0 

7 = h + P (16) 
dy dy 

or f as expected, the applied shear is the sum of the web 
dP 

shear h and the shear carried by the cap strips 

P » . The slope between the cap strips is dh/dy. The 

web shear dP/dy has already been obtained entirely inde— 
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pendently of V "by the methods outlined and equation (ll) 
so that no taper correction will be necessary. 

Equation (16). can be used, however, as a rough check 
on the shears obtained from equation (ll) by conputing the 

tern P -|k fron the axial loads and slope, adding it to 

the vertical web shears, and checking the sun with 7, the 
applied beam shear on the section. 



Tabular Forms, for Shear Computation 

Table 3 has been developed for quick computation of 
shear stresses and shear center. A brief explanation is 
given here, 

For ease in integration, elements should be. numbered 
starting from the upper cap strip in a clockwise direction. 
Shear panels should be numbered in a clockwise direction 
corresponding to the adjacent counterclockwise element. 
Material in flanges extending beyond the torsion box should 
be numbered separately starting from the free edge. (See 
fig. 23.) 

The cut will be taken in the vertical web, web A, 

with thickness t_. Column 2 contains the width ds of 

a 

each shear panel measured between elements. Columns 6 and 
7 are the coordinates x and z of the various elements 
measured to the center of gravity of the section. These 
coordinates differ slightly from those of columns 12 and 13 
of table 1 because in table 1 the coordinates are to the 
center of gravity of each element, while in table 3, the 
coordinates are measured to the point at which the element 
is attached to the skin. This method is necessary in order 
that the summation of moments will be accurate. These coor- 
dinates may be scaled off the rib drawings. 

Columns 8 and 9 should be watched for signs and the 
absolute value of their algebraic difference should be 
entered in column 11. Column 10 need be used only where 
there is sharp curvature in the skin between elements, such 
as in the nose D. Columns 6 to 11 may be omitted if the 
shear-center calculation is neglected. 

Columns 12, 13, and"T4 give' the" HW" of" change of ■ ax- 
ial load for each element at the station in question. A 
tangent to the load curve is drawn at the station (see fig. 
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24) and the slope is computed. Column 15 is the progres- 
sive summation of the values in column 14 starting on the 
clockwise side of the cut vertical web. In column 17, q], 
is the shear stress in any member due to bending only. (See 
equation (ll).) From the shear-center location (calculated 
or estimated) and the computations in lines 3 to 9, the tor- 
sional shear q t is obtained, which is added to q D giv- 
ing column 18. The right-hand term of equation (15) is ob- 
tained from the summation of column 19 and is used in the 
shear-center calculation in lines 12 and 13, 

If it is desired to omit the shear-center calculation 
and solve for the shear stresses by the statically determi- 
nate method of references 7 and 10, the tabular form may 
still be used with slight modification. Columns 2 to 5 may 
be omitted. Columns 16, 24, and 30 should be changed to 

E ^ x 2AA, [column (l5) x column (9)3, Columns 17 and 

19 may be left out and a separate computation may be made 
for the correction shear at the cut section 

7 B x c.g. ~ V C z c.g. - column (16) 

2 x area of torque tube ^b+t 

which is then added to column (15) to get q in column 
(18). 

Table 3 is worked out for station. 232 of the wing for 
airplane B as an example. This section has projections be- 
yond the torque box, which are handled as previously noted. 
Care must be taken that the proper sign is used for these 
projection shears. Note that the top-cover projection 
changes sign when added to the torque box in column (15). 

For simplicity, the shear from element V is assumed 
to be divided evenly between top and bottom cover rather 
than computed as though it were a double torque box. 
Shears in the projections are unaffected by the torque 
shears on the box. 

The load curve of figure 24 is a good check on the 
general trend of the stresses and the loads. Each 
stringer or cap strip should show a smooth curve. The slope 
of each curve at any point represents the load per inch that 
the rivets must transfer from skin to stringer. This fact 
can be used in designing the rivet spacing. The spacing 
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will "be most critical where the curves are steepest , that 
is, at each end of the stringer. 

•The curves for J and G are caused "by the reverse 
bending in the front "beam, and the irregularities through 
the wheel-well region (station 142) are apparent in the 
cap strips B and B". 

The values of AP/Ay of column 14 were obtained by- 
drawing tangents to the load curves at station 232,. reading 
values of P at station 222 and station 242, and dividing 
"by Ay = 20. As a check on the accuracy of the calcula- 
tions, the summation of column 14 should he approximately 



zero . 



ANAL YT I OAT. CHECK FOE EFFECTIVENESS CURVES 



Final proof of the accuracy of the effectiveness curves 
and the axial and shear stresses obtained from them can be 
made by using the principle of consistent deformations. 

In figure 25 assume that the applied loads cause the 
structure to deform as shown. If the effectiveness curves 
and the computed stresses are accurate, then the deforma- 
tion of point H to H' should he the same, whether fig- 
ured "by the path A E» F 1 G 1 H' or hy the path A B 1 C» D 1 
H 1 . The first path will include axial deformation in AE 
plus shear deformations in panels EF, FG , and GH. The 
second path is made up of shear deformations in AB , BC, 
and CD, and axial deformation in DH . 

Axial deformations will be computed hy the formula 

3 = fh TT 
p b e 

where f-^ is the average axial stress in a length Ay of 
a memoer. Shear deformation is computed from the relation 

d = 5 ds 
8 t G e 

The effective shear modulus G Q should be so selected that 
diagonal tension is considered. 



In the following section a check is made for the wing 
of airplane B between stations 23 2 and 277. (See fig. 26.) 
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Calculations for shear stresses at station 277 are not in- 
cluded out are similar to those for station 232. For now 
work, approximate checks, which will serve as a very good 
control, can "be made when the effectiveness curves are 
first drawn. 

TOP COVER 



In order to take diagonal tension into account from 
reference 2, & e will "be selected at 0.3& = 0.2 x 
4,000,000 = 800,000 pounds per square inch. 



SHEAR DEFORMATIONS 



Panel 


Station 232 


Station 277 




q 

(lh/in.) 


t 

(in.) 


ds 
(in.) 


qds 


(Id/ in.) 


t 

(in.) 


ds 
(in.) 


qds 
t&e 


0 
E 

a 

H 


•4.9 
9.8 
19.2 

22.1 


o.oUo 
.oHo 
.oUo 

.OHO 


7. Ho 
11.20 
12.00 

7.00 


-0.00113 

.OO3U3 
.00720 
.00l|S3 


-6.5 
2.0 

19.5 


o.oUo 
.0U0 
.oHo 


15.00 
9.00 
7.00 


-O.OO3OU 
.00060 
.00U27 



Total .OIU33 in. 



.001S3 in. 



AXIAL DEFORMATION 



Ele- 
ment 


Average 
P 
(lh) 


Average 
area 
( sq in. ) 


Ay 

(in.) 


PAy 
AE 


B 


-143 0 


0.404 


45 


-0.01592 


E 


-510 


-.210 


45 


-.01092 


G 


-340 


.204 


45 


-.00750 


H 


-65 


.083 


45 


-.00352 



Then from figure 26, 



Path BB ~ BEG-H 0.01592 + 0.00183 = 0.01775 in. 

Path BCEOH - HH 0.01433 + 0.00352 = 0.01785 in. 

This result checks very well, supporting the validity of 
the effectiveness curves. Such checks can be made wher- 
ever considered necessary. 

Crruraman Aircraft Engineering Corporation, 
Bethpage, N. T. 
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APPENDIX A 
SIGN CONVENTIONS AND SYMBOLS 
Signs 

Applied loads: 

Shear loads - positive up, rearward and outboard 

Bending moments - positive if causing compression in 
upper rear corner of section 

Pitching moments - positive if a stalling moment 

Stresses: 

Axial stresses - positive if tensile 

Shear stresses - positive if in same direction as stall- 
ing moment (clockwise if section is viewed toward 
plane of symmetry of airplane) 

Coordinates: 

x * y» z - positive if rearward, outboard, and up from 
origin 

Elemental skin widths, ds: 

Integrations in clockwise direction (same direction as 
positive shear) 

Symbols 

The symbols used in this paper correspond to similar 
ones in references 1 and 11. 

A cross-sectional area 

Ay is flange, or cap- strip, area and includes any- 
adjacent material considered to act with it. Aj, is 
area of stringers from flange to center of panel and 
includes the effective adjacent sheet. If there is 
""""oUTy'one stringer, A^ is one-half the area of the 
stringer, Aj = Ajp + Aj, 



S3 

P axial load 

0 " shear modulus, including effect of tension fields 

t thickness of cover sheet 

» 

3 Young 1 s modulus 

o half-width of "beam for single-stringer beams; distance 
from flange to center of gravity between flange and 
center line for a mult i stringer beam 



shear-lag parameter J ~ (~ + ~) 



8 shear force 
h depth of beam 
f stress 

2* effectiveness factor (f L /fy) 

S|g and B*q effectiveness factor for beam and chord bending 
M bending moment (Ph) 

I n , I zz t I xs geometric moments of inertia 
V shear strain 

q running shear, pounds per inch of skin 
I total change in length of flange 

L length of panel from free end to plane of symmetry. 
Also total change in length of stringer 

x distance along x axis measured from origin halfway 
between flanges of panel (equation (4~)) t or from 
reference point 

c distance along % axis measured from reference point 

x b c • ys.c.i B s.c. 4 iB*a a c« froa center of gravity of 
material to shear center 

x„ - . y„ „ i s* . distance from load axis to e.g. of 
section 
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V auxiliary parameter 
a angle of attack 

V total work 
I torque 

V applied shear on the section 
V X ,V 2 applied Dean chord shear on section 
as v^idth of shear panel 

9 p axial deformation 

a shear deformation 



Subscripts 


t 


clue t o uoxo[u.c 


b 


due to "bending 


a 


at section a 


e 


effective 


3 


beam 


c 


chord 


F 


flange 


L ' 


stringer 


s 


shear 


3! 


total 


V 


web 


0 


cover sheet 


e.g. 


center of grav 


s . c . 


shear" center 


n 


at section n 


0 


origin 
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APPENDIX B 
DITAILSD DIEIVAIIOH 07 
SPAS VISE EPF35CTIVEHSS8 EQUATION 

The fundamental equations are: 

dP L = dSj, (la) 

dPj c !vj£ . dS 0 (ID) 
n 

af b = if (*i - *L> 4y (Ic) 

In equation (la) 

dP L = df^ dS 0 = f,tdy (a) 

so that if equation (la) 1b differentiated and equation 
(lo) is substituted, 



f- = tt — = - _ kiip - *l ' 

dy dy Al S1jA l 



(H) 



Prom equations (la) and (a) 

<•> 

Iron equations (lo) and (a) 
and from equations (c) and (d) 
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Substitute equation" (e) in equation (b) and simplify 

a 3 f L _ Gt (_1_ j\ St p Sydy 

r s Eb \A L Ap/ L EbA L J A f h 



dy" 



Since 



and 



Let 



Gt / 1 1 \ Gt 1 A L + A F r S w dy 

Eb W + Aw/ f L " Eb AtAtp At. + A™ . / ~~h~~ 



AjAp Ajj Ap 

Al + Ap = Aj 



dx Eb \A. L Ap' n J A T h ' 



Eb \A L Ap- 
Also S'jfdy = dM, so that equation (f) becomes 

il£j - K 3 f L + k 3 JL = o ( g ) 

dy 3 11 A T h 



and by analogous derivations, 



- K ^ * K * iSe - 0 (h) 

£*JL _ K 3 f s + K s !^ = 0 ■ (i) 

dys htA T 



In the NACA methods of analysis, equations (g), (h) . 
and (i) are solved by means of arithmetic integration so 
that if K, M, S w , h, t , Ai, , Ap , and b are variable along 
the span, they may be taken into consideration. (See ref- 
erences- 1- to 6.) - 

The solution of equations ( g) , (h) , and (i) will be 
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obtained for the following conditions: 

E s = a constant 
= a constant D 

S V^L _ a constant R 
htAj 

Let p denote the operation d/dy; then equation (g) be- 
come s 

<p s - **> *i - - sfc 

Take the first derivative to eliminate the constant 

p (p 3 - X s ) f L = 0 (j) 

This is a linear differential equation that has a so- 
lution in the form 

m x y m ? Y m^y 

f L ss Ae . + Be + Ce (k) 

where , m s , and m 3 are the roots of equation (j) for m 
substituted for p and fr = unity. 

m(m 3 - K 3 ) = 0 

or 



or equation (k) becomes 



f L = Ae~ Ky + Be Ky +0 (l) 

If equation (l) is substituted in equation (g) and is then 
solved for 0, it will be seen that 

o = JL = D 

hAj 

Thus the solution of equation (g) for the condition stated 



1 s 



f L = Ae"* 7 ♦ Be Ky + D (m^ 
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Similarly 

fp = Ee- Ky + Pe Ky + D (m 3 ) 

f g = Ne-^ + 0,^ + R U 3 ) 

Coefficients A, B, B, P, IT , and ft are determined by 
the boundary conditions f or . the particular problem con- 
sidered. 

Stringer Btress .- For the beam of figure 3 the fol- 
lowing boundary conditions hold when 

y = 0 f L = 0 

and when ■ df T 

y = L _J± = 0 
dy 

Substitute in equation 

A + B + D = 0 

-KAe~ KL + KBe KL = 0 

Solve simultaneously for A and B 

A = " De 



2 cosh KL 



2 cosh KL 
Substitute in (m x ) and simplify 



U = » 



|~ K(L-y) -K(L-y)J 



L ~ 2 cosh KL 

f L = D fi - cosh g(I-y) T (2b) 
{_ cosh KL J 

OaP-strixi stress .- The boundary conditions for the 
cap strip are : 
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"When 

y = 0 



and when 

y = L 

= 0 



dy 

In a manner similar to the stringers, 

■p - Ti [~i 4. Al cosh K(L~y ) ") / % 

f T '= D L 1 + — cosh EC — J (2aJ 

APPENDIX C 
DEEIVATION OF EQUATIONS FOE BEAM SHEAS 
IN DOUBLE -CEIL TOEQUE BOXES 

In terms of the shears (see fig. 27) at the cut sec- 
tions q_x and q_ s , the variable shears in the three 

parts of the loop q-^ , q t(J , and are 

dP b dP 

q h = *i ~ I d? f s^ = f s x *i - f dy 

i , s ' c 1 3 i f s dy 

d d P d 
a , = <1 + 2 f_ t = f f . t a + 2 — 

M i 3 s dy s d & 3 3 g dy 



In order to find the values of and q g , the 

least-work solution is applied (see reference 12) to the 
two loops (l) and (2): 

c a 
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15 C 

V 3 = ^ Jl f .» td. + ^ f f Sc 2 tds (b 3 ) 



a 

Perform the operations 



9W_ 3W_ 



and equations (t^) and (b s ) become 



c 1 a 



But 



df B " t * df_ = t • bf a = t ' oT7~ = t 
s i s i s s s s 

Substitute these expressions in equations (c x ) and (c 3 ) and 
simplify 



b f s lj t i ds + J f s t i ds = 0 ^ d i^ 
c 0 a ° 



d f« t s ds - / f, t 3 ds = 0 (d-J 



8,2" / S "2 

d «_/ c 

a 



Substitute equations (a) in equations (d x ) and ( d 3 ) 



.G .... .., .. „ 

t x c dp 

EL 
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q 2 ds 



Simplify 




- / — ds 2 ^ (e x ) 





d dP 

£ £ ' <e a ) 



2 



dy 



Equations (ej.) and (e 2 ), although apparently compli- 
cated, are easily obtainable in a tabular computation by 
progressive summation. Values of q_j. and q s can be 
solved for simultaneously and substituted in equations 
(a) for the shears at any section. 

c dp 

In the integration i E g for the common web, both 

upper and lower cap strips should be involved, while in 
the other integration they should be omitted. Skin effec- 
tive in bending may be included in an approximate manner 
by considering the area concentrated at the center of the 
panel and then treating it as an extra stringer. 

Triple torque boxes may also be handled by a similar 
derivation, there being three equations such as equations 
(e 1 ) and (e s ) to solve simultaneously. 



The shear center may be found in the same manner as 
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for a single torque box as soon as the bending shears are 
determined, by taking a moment of these shears about the 
center of gravity and equating it to the moment of the ex- 
ternal shears about the center of gravity. 
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Figure 12.- Effectiveness factors for box beam for airplane-A wing. Center sectionjtop surface, to 
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HACA method is for beam bending only. Chord bending would raise front-beam stress and lower rear-beam 
stress without changing curve shape and would improve agreement. 
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Figure 13,- Chordwise 1 g stress distribution. Station 43, Airplane— A wing "box beam;top surface; a=20 
NACA method is for beam bending only. Chord bending would raise front-beam stress and lower rear-beam 
stress without changing curve shape and would improve agreement. 
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